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Abstract. A nonlinear Schrodinger equation arising from light propagation down 
an inhomogeneous medium is considered. The inhomogeneity is reflected through 
a non-uniform coefficient of the non-linear term in the equation. In particular, a 
combination of self-focusing and self-defocusing nonlinearity, with the self-defocusing 
region localized in a finite interval, is investigated. Using numerical computations, the 
extension of linear eigenmodes of the corresponding linearized system into nonlinear 
states is established, particularly nonlinear continuations of the fundamental state and 
the first excited state. The (in)stability of the states is also numerically calculated, 
from which it is obtained that symmetric nonlinear solutions become unstable beyond 
a critical threshold norm. Instability of the symmetric states is then investigated 
analytically through the application of a topological argument. Determination of 
instability of positive symmetric states is reduced to simple geometric properties of the 
composite phase plane orbit of the standing wave. Further the topological argument 
is applied to higher excited states and instability is again reduced to straightforward 
geometric calculations. For a relatively high norm, it is observed that asymmetric 
states bifurcate from the symmetric ones. The stability and instability of asymmetric 
states is also considered. 
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1. Introduction 

Inhomogeneities can act as an effective trapping to collective excitations in nonlinear 
media. In the field of nonlinear integrated optics, the first theoretical works on guided 
waves by an interface between a linear and nonlinear medium appeared in [T71 dj [29j [18] . 
Various stationary wave profiles propagating along nonlinear planar optical guides in 
a layered structure are then extensively considered (see, e.g., [2j [TlJ [23j IU [26] and 
references therein). When the inhomogeneities are periodic, one will obtain discrete 
waveguide arrays, which have become an independent topic of interest |15j . A next 
fundamental question is whether or not the standing waves are stable to propagation 
along the nonlinear waveguide. 

The stability of stationary nonlinear Schrodinger waves in homogeneous media was 
first considered by Vakhitov and Kolokolov [3H [H]. Using variational arguments, 
a criterion was derived relating the soliton linear stability and the slope of the 
corresponding power-dispersion curve, i.e. the known Vakhitov-Kolokolov condition. 
The method was later rigorously justified by Weinstein in [35], and in [9]. Subsequent 
studies extend the condition for various situations, including inhomogeneous problems 
[T2l 1301 131T f32l 1331 [19] (see also a recent brief review [25] and references therein) . 

An interesting waveguide system was proposed in [31], consisting of a self-focusing 
Schrodinger equation and a self-defocusing type inhomogeneity with finite length. It 
is experimentally feasible to fabricate such a waveguide using the current technology 
as self-focusing and self-defocusing can be achieved in the same medium, structure, 
and wavelength [20]. In the context of a Bose-Einstein condensation [61 [3J, [5], 
which is also modeled by a nonlinear Schrodinger equation [10J [2TJ, such a sign- 
changing nonlinearity coefficient can be created by spatially varying the condensate's 
atomic scattering length making the so-called collisionally inhomogeneous nonlinearity 
[271 EH [22]. 

Tran [31] extended the work of, e.g., [21 [12], in which the inhomogeneity is linear. In 
[3Tj the self-defocusing inhomogeneity has a small nonlinearity coefficient, such that the 
characteristics of the stationary solutions are closely related to the corresponding linear 
problem. The system is later studied by Leon [16] . where the nonlinearity coefficient of 
the inhomogeneity is of the same order as the self-focusing regions. Analytical solutions 
of stable symmetric solutions below a threshold amplitude are derived in terms of 
Jacobian elliptic functions [16]. Here, we revisit the problem. 

We study the existence and stability of symmetric and asymmetric solutions, 
particularly the fundamental and the first excited mode, when the nonlinearity 
coefficient of the defocusing inhomogeneity is of the same order as the focusing bounding 
regions. We show that continuing from the linear limit solutions, there is a saddle-node 
bifurcation at which the symmetric mode becomes unstable and asymmetric modes 
emerge. Even though it is similar to the results reported in [31], there is a significant 
difference where the asymmetric positive solutions are all stable in their existence region. 
Moreover, [31] only considers positive solutions. Besides determining the instability of 
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symmetric solutions numerically, we also show it analytically using topological argument 
techniques as developed in [12j 02] ■ We also comment on the inapplicability of the 
analytical techniques to asymmetric solutions. 

In Section 2, the governing equations are discussed and the corresponding linear 
eigenvalue problem is derived. In Section 3, we consider the linear limit of the equations, 
where a transcendental equation determining the bifurcation points of nonuniform 
solutions from the uniform solution u = is derived. In the same section, numerical 
continuations of the fundamental and the first excited state from the linear limit to 
nonlinear states are presented. The linear (in) stability of the numerically obtained 
(symmetric and asymmetric) solutions is then determined numerically by solving the 
corresponding linear eigenvalue problem. The instability of the symmetric solutions are 
analyzed analytically in Section 4 using a topological argument. In Section 5 we consider 
some asymmetric solutions. 



2. Mathematical model 



We consider the following governing system of differential equations 

i% + *,x + |^| 2 ^ = \x\>L, 
i^t + - r]\fy\ 2 fy = \x\ < L, 

where the 'outer' and the 'inner' equation has focusing and defocusing (77 > 0) type 
nonlinearity, respectively, and L is a positive real parameter representing half the length 
of the waveguide. The norm 

/oo 
\^(x,t)\ 2 dx, 
-oo 

which is physically related to the intensity power of the electromagnetic field in the 
context of nonlinear optics or the number of atoms in Bose-Einstein condensates is 
conserved. 

To study standing waves of (JTJ), we pass to a rotating frame and consider solutions 
of the form ty(x,t) = e"~ tult ip(x,t). We then have 

+ i>xx + \i>\ 2 ^ = (V - u)ip \x\ > L, 
i*l>t + ipxx - ?7 M 2 ^ = \x\ < L. 

Standing wave solutions of flTJ will be steady-state solutions to (j2J). In the following, 
the parameter 7] is taken to be rj = 1. We consider real, t independent solutions u(x) to 
the ODE: 

u xx = {V-u)u-u 3 \x\>L, 
Uxx — —uu + u 3 \x\ < L. 

To obtain solutions that decay to as x — > ±00, the condition that V — u > is 
required, with u G R. We will also require that u x — > as x — > ±00. To establish the 
instability of a standing wave solution we linearize (j2J) about a solution to (j3J). Writing 
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ip = u(x) + e ((r(x) + is(x))e xt + (r(x)* + is(x)*)e A **) and retaining terms linear in e 
leads to the eigenvalue problem 

, / \ / D_ 

A 



r 



M 



s / \ -£>+ 
where the linear operators D + and Z?_ are defined as 

9L-(V-u) + 3u 2 , \x\>L, 



r 
s 



a 2 

dx 2 

9a; 2 
jS 2 



u) — 3u 2 , 



\x\ < L, 



+0J-U 



(V — id) + M 2 , |x| > L, 

2 Id < L. 



(4) 



(5) 



(6) 



It is then clear that the presence of an eigenvalue of M with positive real part implies 
instability. 

3. Linear states and their continuation 




Figure 1. Bifurcation points of non- uniform solutions from the zero solution u = in 
the (L,w)-plane for V = 1. The insets present a sketch of the corresponding solution 
u(x) along the first two branches. 



In the small limit of u(x), the governing equation (jHJ) is reduced to the linearized 

system 

Uxx = (V — u)u \x\ > L, 
u xx — —uu \x\ < L, 
which can be simply solved analytically to yield 



(7) 



-Vv-m\x\ 



x < —L, 



u(x) 



c e cos(y/ux) + c Q sm(y/uix), \x\ < L, 

x > L. 



(8) 



From the natural continuity conditions at the points of discontinuity 
u(±L+) = u{±L-), u x {±L+) = u x {±L~), 
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one will obtain that the parameters of the linear states above will have to satisfy the 
transcendental equation 




This equation determines bifurcation points of non-uniform states from the zero solution. 
A plot of ([9]) for V = 1 is given in Fig. [TJ 

Starting from a bifurcation point, as the parameter u varies, the corresponding 
linear limit solution will deform and nonlinear terms will play a role. Even though one 
can still represent the continued solutions in terms of the Jacobian elliptic functions 
|16j . here we solely use numerical computations. A pseudo-arclength method is used 
to follow the existence curve of a solution as a parameter is varied. We have solved 
Eqs. (J3j) — (HD numerically to study the existence and the stability of localized standing 
waves, where a central finite difference is used to approximate the Laplacian with a 
relatively fine discretization. In particular, we consider the first two branches of linear 
limits shown in Fig. [TJ 

3.1. Positive solutions 
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Figure 2. (Color online) The norm P as a function of uj for the first (positive) state 
corresponding to the first branch in Fig. [TJ As u; increases from the bifurcation point 
uj R3 0.265, there is a bifurcation at which the state becomes unstable. In addition to 
the symmetric positive states, there is also a stable asymmetric state along the middle 
(solid blue) branch. Solid and dashed curve represents stable and unstable solutions, 
respectively. 

First, we consider the continuation of the linear state corresponding to the first 
branch. For illustrative purposes, we take V = 1 and L = 2, i.e. the fundamental 
state mode originates from u ~ 0.265. In Fig. [2j we present the numerically obtained 
continuation of the linear positive solution as u varies. 

At the bifurcation point, the linear state is expected to be stable, similar to the zero 
uniform state u(x) = 0. As w increases, the norm iV of the solution increases as well. 
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By appealing to the work of [HI I2H E2] , we obtain that the solution along this branch 
is stable. In Fig. Eta), we depict a solution corresponding to point A in Fig. [2] and its 
eigenvalue structure in the complex plane, where one can see that all the eigenvalues 
are on the imaginary line. 
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Figure 3. Profile of solutions at points indicated as A-D in Fig. [2] and their eigenvalue 
structures in the complex plane. The phase-portraits of the solutions in each panel are 
presented in Fig. IA11 



As the parameter u is increased further, there is a bifurcation at which the existence 
curve reverses direction. Symmetric solutions are unstable along this branch. Two 
solutions and their eigenvalues in the complex plane corresponding to point B and C 
are shown in Figs. Etb) and (c), respectively. One can note that the instability of the 
solutions are due to the presence of a pair of eigenvalues with nonzero real part. 

Interestingly, in addition to the symmetric states, at the bifurcation point where 
symmetric states become unstable, there is an existence curve emerging, corresponding 
to some numerically stable asymmetric states. A solution indicated as point D in Fig. 
[2] is presented in Fig. [3](d) . 

Figure [2] is similar to Fig. 1 in (31], which is for the case of relatively small rj > 0. 
Writing (V — ui) —> W 2 and —oj — > {W 2 — V 2 ) and L = 1, the time independent system 
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of (j2J) becomes the same as Eq. (1) in |31j. One important difference is that in our case, 
asymmetric solutions are all (at least numerically) stable in their existence domain. 

When a solution is unstable, it is certainly of interest to see the dynamics near it. 
Here, we have solved the time-dependent governing equation (JTJ using a Runge-Kutta 
method. Depicted in Figs. H^a) and (b) are the dynamics of solution (b) and (c) in Fig. 
El respectively, perturbed initially by small random disturbances. Shown is the modulus 
\i/j(x, t)\ 2 . One can clearly see that the instability of solution in Fig. [3^b) manifests in 
the form of spontaneous symmetry breaking, while the instability of the solution in Fig. 
[3](c) is in the form of a soliton generation, similar to that reported before in |16j . 

3.2. First excited state 

We have considered as well solutions bifurcating from the first excited corresponding to 
the second branch in Fig. [TJ For the same V and L as above, this state bifurcates from 
the point u ~ 0.898. We depict in Fig. |5]the continuation of this state as to varies. 

As the parameter u increases from the bifurcation point, one will obtain a 
numerically stable symmetric state. When the parameter is increased further, there 
will also be a 'direction reversal' point, where the symmetric state becomes unstable, 
similar to the case of fundamental state solutions above. Shown in Fig. [6](a) is an 
example of this state and its spectrum in the complex plane. 

Besides similarities with the previous case, we also observed several differences 
here. These include the fact that there are now more than one existence branches 
corresponding to asymmetric solutions. Presented in Fig. E^b) is a solution along the 
first asymmetric branch, indicated as point B in Fig. [5j and its eigenvalue structure. As 
a comparison, we also plot in Fig. [6][c) and (d) the symmetric and asymmetric solution 
from point C and D in Fig. [51 respectively, and their eigenvalues in the complex plane. 
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Figure 5. The same as Fig. [21 but for the first excited state, corresponding to the 
second branch in Fig. [TJ The bifurcation point of the linear state is uj w 0.898. 

The two asymmetric solutions from point B and D above are clearly different. 
By viewing this first excited state as composed of two static out-of-phase solitons, the 
asymmetric solution B can be seen as composed of two solitons with one of them spatially 
displaced, while the solution D can be viewed as composed of two solitons with different 
amplitude. 

It is also interesting to note that the asymmetric solution D is not always 
numerically unstable in its existence region. Numerically, we observed a region of 
stability and of instability for this asymmetric state, depicted as a solid and dashed 
line respectively in Fig. [5j In Fig. Ele), we present an unstable asymmetric state and 
its spectra, where one can see the presence of two pairs of eigenvalues with nonzero real 
part. 

Typical time dynamics of unstable solutions for this case is presented in Fig. [71 
In particular, we plot the dynamics of asymmetric solution B and E under random 
perturbations. The dynamics in time of solution C is similar to that of solution B. 

4. Unstable symmetric solutions past the bifurcation point 

In the following, we will analytically prove the instability of symmetric solutions past 
the 'direction reversal' point above. To show instability of the standing waves, we will 
show that M has a real positive eigenvalue. This is done by applying the main theorem 
of [IB]- It can be shown that the following quantities are well defined (see for example 

my- 

P = the number of positive eigenvalues of D + 
Q = the number of positive eigenvalues of D- . 
We then have the following: 

Theorem 1 ([H]). If P — Q ^ 0,1, there is a real positive eigenvalue of the operator 
M. 
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Figure 6. Profile of solutions and their corresponding eigenvalue structure in the 
complex plane at points indicated as A-D in Fig. [5] The solution phase-portraits are 
shown in Fig. IA2I 
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From Sturm-Liouville theory, P and Q can be determined by considering solutions 
of D + v = and D_v = 0, respectively. In fact, they are the number of zeros of the 
associated solution v. Notice that D_v = is actually satisfied by the standing wave 
itself, and that D + v = is the equation of variations of the standing wave equation. It 
follows that: 

Q = the number of zeros of the standing wave u(x). 

P = the number of zeros of a solution to the variational equation along u{x). 

We will focus first on the case when we have a positive localized steady state solution, 
i.e. Q = 0. 



4-1. Positive solutions 

The idea is to use a dynamical systems point of view, and geometric properties of the 
solution curves in the phase portrait to establish when P > 2. The t independent 
solutions to equation <^ can be represented by composite phase portraits constructed 
by a superpositioning of the phase portraits of the 'outer' system: 

Ux = y, Vx = (V - u)u - u 3 , (10) 

and the 'inner' one: 

u x = y, y x = -uu + u 3 . (11) 

We can view the composite picture as a single, non-autonomous system with phase plane 
given by: 

Ux = y, 

(V-u)u-u 3 , \x\ > L, (12) 
— oju + u 3 , \x\ < L. 
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In the phase plane of (1101) . the outer system admits a soliton solution, given by the 
equation: 

y 2 =(V-u)u*-^, (13) 
while the inner system ( TTTj) admits a heteroclinic orbit in the phase plane given by: 

2 2 u4 u2 /-.as 

y 2 = -uu 2 -- + —. (14) 

The conditions that u and u x decay to zero as x goes to ±oo, mean that in the superposed 
phase portraits, all steady-state solutions that we are interested in will lie on the soliton 
curve as x — »■ ±oo. 

Solution curves of the inner system are given by: 

y 2 = -cou 2 + — + a (is) 

The solutions we are interested in will travel in the phase plane out from the origin along 
the homoclinic orbit of the outer system described by ( TT3l) and then 'flip' to the inner 
system for 2L units of x, and then 'flip' back to the outer system along the homoclinic 
orbit. Define (uq, yo) as the point in the phase plane of (f!2l) where the solution initially 
flips from the outer to the inner system, and define as the point in the phase 

plane where the solution returns to the outer system. Using this notation we can write 
the equation of the part of the solution curve in the inner system in the (composite) 
phase plane as: 

y 2 = -uju 2 + ^ + Vu 2 - 4. (16) 

As the parameters, V and u vary, the relative position of the homoclinic orbit 
described by ([13]) and the heteroclinic orbit described in (TT4"]) will change. We consider 
the qualitative and numerical differences that occur for various values of L corresponding 
to different configurations of the potential V and the frequency u. Namely when the 
curves described by (JHJ) and ( fT3|) are 'close' together, linearly unstable standing wave 
solutions to JT]) appear. 

Theorem 2. Unstable positive localized solutions to (T7]j occur whenever 77 < §■ The 
unstable symmetric solutions are the ones which are symmetric in the phase plane 
with respect to the u-axis and leave the homoclinic orbit at a point (u ,y ) satisfying 

The theorem will be proved by showing that P > 2 for such a solution. From 
Sturm-Liouville theory we have that P will be the number of zeros of a solution to the 
variational equation D + v = satisfying the boundary conditions v — > as x — > — oo. 
A solution to the variational equation associated with a solution of ([3]) can be found by 
following a tangent vector around the orbit under the flow of the linear equations ( Tl2l) . 
Satisfying the initial conditions means the solution will be a tangent vector to the orbit 
of the solution of ([3]) in the phase plane until the discontinuity of (J3]). The number of 
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zeros of such a solution can be found by determining the number of times that such a 
vector must pass through the vertical as the base point ranges over the entire orbit. 

Denote by b°(u, y) = (b 1 (u, y), b 2 (u, y)) a tangent vector to the homoclinic orbit 
of the outer system at the point (u,y), and let h\u,y) = (b{(w, y), h^u, y)), be a 
tangent vector to an orbit of the inner system at the point (u, y). As (u, y) varies along 
the homoclinic orbit, h°(u,y) = (p, q) solves the linear system: 



(17) 



( 18 ) 



Px = q, 

q x = (V - tu)p - 3u 2 p 
while b 7 (it, y) solves the linear system: 
Px = q, 

q x = —up + 3u 2 p 

as (u, y) travels along a curve given by ffl5|) The composite linear variational equation 
is then given by 

Px = q, 

(V-u)p-3u 2 p, \x\ > L, (19) 
— up + 3u 2 p, \x\ < L. 



We denote the explicit vectors, h°(uo, yo) and b J (ui, yi), by j3° and /3 1 respectively. Let 
F° and F 1 denote the flow of the outer (flOj) and inner (TTTj) systems respectively and 
$° and $ 7 denote the variational flow of the outer and inner systems along solutions. 
If d(u, y) = (di(u, y), d 2 (u, y)) is a vector in the tangent space of the phase plane at the 
point (u,y), and if (u, y) is flowed under F° and F 1 to F°(u, y) and F T (u, y) respectively, 
let $°(d)(F°(u, y)), $ 7 (d)(F / (w, y)) denote the image of the vector d under the flow 
of the outer and inner variational systems respectively. 

Before moving to the proof of the theorem we need two geometric facts about the 
flows $° and $ 7 . 

Fact 1. The variational flows $° and $ 7 map the tangent line to a solution to / TiTjj) 
and ( 07]) at a point (u, y), to the tangent space of the solution at the point F°(u, y) and 
F^u^y) respectively. 

Fact 2. The flows $° and are orientation preserving. 

This second fact is reflected in the following way; if b(u, y) = (bi(u, y), b<i{u, y)) and 
d(u,y) = (di(u,y),d 2 (u,y)) are two vectors tangent to the phase space at (u,y), then 
the sign of the cross product of the two vectors is unchanged under the flows. That is: 

sgn(h(u,y) x d(u,y)) = sgn($°(b)(F°K y)) x ® (d)(F°(u,y))) (20) 

For every point in the orbit of the flow F . The above also holds with the outer flow 
replaced by the inner flow. 

A consequence of these two facts for the system of interest is the following lemma: 
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Lemma 3. If d(u, y) is a vector in the tangent space to the phase plane of the homoclinic 
orbit at the point (u,y), with u > such that d(u,y) x b (u,y) > 0, then 

Km $°(4)(J^(«, *)) = *( J— ) (21) 
i a\ \ y V — co I 

where k is a positive real number. 

Proof. Because b (u, y) is a tangent vector to the homoclinic orbit, and u > 0, 

Um<& (^)(F°(«,y))= hm^-(F°(u,y)) ' / ' 



ib°r v " *-+oo|b°p v 11 vi + v-uj\ vv-cu 

and since d is not tangent to the homoclinic orbit, 

lim $°(-±-)(F°(u,y)) = k( ±- |, (23) 
x^oo |d| \ v V 7 — co 

where k is some real number. This is because ^ ) is the unstable asymptotic 

eigenvector of (TIT]) . But from the first and second fact we have that 



< lim * (±)(F°(u,y)) x ±(F°(u,y)) = ^^L , (24) 

X^-OO |d| |D| y/l + \/ — 

which means that must be positive. □ 

In order to determine the number of zeroes of a solution to D + v = 0, and v — > as 
x — > — oo, we need to determine the number of times a vector gets pushed through the 
vertical by the flow as its base point moves along the orbit. In order to count this, we 
break the orbit up into three parts. These are when x < —L, \x\ < L and x > L. Let A\ 
denote the number of zeros of the solution to the variational equation as the basepoint 
corresponds to the range x < —L, A 2 the number for \x\ < L, and A 3 for x > L. We 
have the following: 

• The number A\ is the number of times b (it, y), a vector tangent to the homoclinic 
orbit passes through the vertical as (u, y) travels from (0,0) along the homoclinic 
orbit of the outer system to the point (ito,2/o)- 

• The number A 2 is the number of times $ 7 (/3°)(w, y) passes through the vertical as 
(u, y) travels from (w , Ho) along the path y 2 = ^ — cou 2 + Vu\ — u\ to the point 
Oi,?/i). 

• The number A 3 is the number of times $ ($ 7 (/3 yi))(u, y) passes through the 
vertical as (u, y) travels along the homoclinic orbit of the outer system from (ui, y\) 
to (0,0). 

The number of zeros of a solution of (fT9]) will then be P = A\ + A 2 + A 3 . Theorem 
[Uthen says that if u(x) > and if A\ + A 2 + A% > 2, then the underlying orbit in the 
phase plane represents an unstable standing wave solution of ([I]) . We are now ready to 
prove Theorem [2] 
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Figure 8. (a) An unstable orbit of the first type. The dashed line is the homoclinic 
orbit of the outer system. The dot represents the jumping off point. The black arrows 
represent the image of a tangent vector to a solution under the variational flows Q ' 1 
at various values, (b) An unstable orbit of the second type. The dot represents the 
jumping off point. The black arrow represents where $(b°) stops being tangent to 
the orbit of the solution. 



Proof of Theorem^ We have two cases to consider, the first is when \j \ < uq < 
and the second is when u Q > y/u. In both cases however uq > \f\T 



ijj. 



Case 1. Here — < w < y/u. An example is as in figure 8(a) . Note that in this case 
Uo > 0. We show that A 2 > 1 and A 2 + A 3 > 2, so the solution represents an unstable 
standing wave. First we show that A 2 > 1. Notice that (3° = (y , (V — uj)u — Uq), 
and that the first coordinate is positive while the second coordinate of (3° is negative 
since ^ < uq < yfu. Now notice that b 7 (w , y ) = (y , — uu ), and again the second 
coordinate u > \/V — u. This means that f3° x b 7 (-u , y ) = y u (2ul — V > 0, because 



u ° > V 2 an< ^ ^° ^ S P os ^^ ve - Now we flow both vectors along the inner system until 
we get to the point where the solution to the original equation crosses the u axis in the 
phase plane of the system described by f fT2|) . Denote this point by (wmax,0). We then 
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have the following: 

< $ / (/3°)(«max,0) x b 7 ( Mmax ,0) = « ax - ^ max )$((/3°)(w max , 0). 



Since u < \fw this means that Umax = y w — \J — 2(y«g + Uq) < y^u; and so the 
first coordinate of $ 7 (/5°)(w m ax, 0), $((/3°)(w m ax, 0) < 0. But this means that the flow 
has pushed the original tangent vector to the homoclinic orbit through the vertical at 
least once by this point, since the sign of the first coordinate has changed, so A 2 > 1. 

To study A 3 , we must determine the number of times that yi))(u, y) 

passes through the vertical as the base point travels along the homoclinic orbit from 
(ui, yi) to (0, 0). We first apply the result of the lemma to the vector ft 1 . Observe that 
(3 1 = (—yo,ul — (jju ), (since in this case (u\,yi) = (u ,—y )) and that the first, and 
second coordinates of j3 are neg ative. Now (3 1 x b°(ui,yi) = yi(Vui - 2u\) > 0, since 

U\ = uq > So applying the result of lemma [3] means that there is a point on the 

homoclinic orbit (u',y') say, where $ (u',y') is pointing vertically upward. Write 
y') := (0,a), with a positive. Now we apply the fact that the variational 
flows are orientation preserving. We have 

< $ ($ / (/3 )K, 2/1 ))K,2/') x (0,o) = a($?($'G8°)(«i, yx))(u',y')) 

which means that the first coordinate of $° ((3°){ui, yi))(u', y') must be positive at 
this point. The first coordinate under the variational flow was shown to be 

negative above, so this means that the vector must have either passed through 

the vertical once more, or been pushed through the vertical by $°, by this point, so 
A 2 + A 3 > 2. 

Thus we have that the number of zeros to the variational equation must be greater 
than or equal to two, so the corresponding orbit must represent an unstable standing 
wave. 

Case 2 In this instance Uq > y/oo. An example of the phase portrait of an orbit of 
this type is in figure 8(b) . Here we remark that since uq > \fu this means that y$ must 
be negative. In this case we show that A\ > 1 and that A\ + A 2 + A% > 2. 

To see that A\ > 1, we note that the tangent line to the homoclinic orbit is 
vertical when the homoclinic orbit crosses the u axis. As the base point moves along 
the homoclinic orbit from (0,0) to (uo,y ), it must cross the u axis, so A 1 > 1. 

Now we study A3. Call the point where the homoclinic orbit crosses the ti-axis 
(amaxi 0). Note that j3° = (y , (V — w)u — Mq) and that the sign of the first coordinate 
is negative. The sign of $ / (/5°)(mi, yi) x h°{ui,yi) can be either positive, negative or 
zero. If it is positive apply lemma [3] to see that 

r $ ($ J (/7 )(m,7/i))(u,y) / k 

where k is positive, so the flow must have pushed the vector through the vertical once 
more, since the first coordinate in the limiting vector is positive. Thus A\ + A 2 + A3 > 2. 
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If $ 7 (/3°)(-Ui, yi) x b°(ui,yi) = 0, then $ / ( ( 5°)(-u 1 , y x ) is tangent to the homoclinic 
orbit at (ui, 2/1). But the flow $ maps tangent vectors to tangent vectors, so 
$ ($ / (/3 )(mi, yi))(wmaX) 0) is vertical in this case, so in this case one can see directly 
that A 3 > 1. 

Lastly if $ / (/3 )(«i, yi) x b yi) < 0, then denote b (-u ma x, 0) := (0, — a) where 
a is some positive number. We then have 

> <f o (<f / (/3°)(« 1 ,2/ 1 ))(«max,0) X (0, -a) = (-a)$f ($ 7 (/3°)K, yi))(u max , 

This means that the first coordinate of the image of the vector under the flow is positive, 
so the vector must have passed through the vertical at some point. Thus AL+A2-I-A3 > 2 
in this case as well. 

So we have that for solutions of this type the number of zeros of the associated 
variational equation is greater than or equal to two. This completes the proof of the 
theorem. □ 



3 




Figure 9. The spectrum > of D+ for unstable symmetric solutions along the 
upper branch in Fig. [5] as a function of w. The inset zooms in on a small region clearly 
showing that there are two curves. 



Remark 1. Solutions shown in Fig. E^b) and (c) (see also the corresponding phase- 
portraits in Fig. IAl( b) and (c)) belong respectively to case 1 and 2 in the proof of 
Theorem [2] above. Let q 6 K. be the real part of a A in the spectrum of D + such that 
D + v(x) = Xv(x). In Fig. [9] we plot the numerically obtained spectral parameter q > 
of symmetric solutions along the upper branch of Fig. [2] as a function of u, from which 
one can see that P = 2 indeed. 

So far we have been primarily concerned with positive solutions to ()3]), but Theorem 
CD applies equally well to standing waves u(x) which are not strictly positive, so long 
as they are smooth enough and both u and u x tend to ±oo. Likewise, the 

techniques used to calculate P in the proof of Theorem [2] do not require that the solution 
be positive. As such, we can apply Theorem (TJ and the techniques above to establish 
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(b) 

Figure 10. (a) An unstable orbit where P = 3 and Q = 1. (b) An unstable orbit 
where P = 2k + 3 and Q = 2k + 1 and k > 1. The number fc describes the number of 
periods the solution stays in the inner periodic orbit. 

the linear instability of some excited states with Q > 1. The first excited standing wave 
we will deal with has Q — 1. 

Fzrsi excited unstable states 

We will consider steady state solutions with the following properties 

(i) The switch from the outer to the inner system takes place in the phase plane at 

v 



(u , y ), with y < and < u < y ~. 

(ii) The solution returns to the outer homoclinic orbit at the point (ui, yi) = (—u , y ). 

Such a solution will be symmetric with respect to the y axis in the phase plane, and 
such solutions exist for all values of V and uj. An example of the phase portrait of such 



a solution can bee seen in figure 10 (a 



Corollary 4. Let u(x) be as above, then P > 3 and so we have a linearly unstable 
standing wave to equation (QJ]. 

In order to prove the corollary, we need the following lemma, which is the analogue 
to lemma [3] in the half of the phase plane where u < 0. 
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Lemma 5. If d(u, y) is a vector in the tangent space to the phase plane of the homoclinic 
orbit at the point (u,y), with u < such that d(u,y) x b (u,y) > 0, then 

J™*°0< FO <^>> = *(^) (25) 

where k is a negative real number. 

The proof is exactly the same as in the proof of lemma [3j except that: 
b°h,„, b° 



hm $°( >)(F°(u,y)) = lim -^(F°(u,y)) 

1 / -J 



-y/l + V-U \ W 



(26) 



because we are in the left half plane. This change of sign therefore changes the sign of 
k. 

It is straightforward to calculate that A\ > 1. This is done in exactly the same 
way as in the proof of the second case of theorem [2j Similarly, as in the second case 
of theorem El one shows that A\ + A2 + A3 = P > 3. This is done by using the cross 
product of the appropriate vectors under the flows of the outer system, showing that the 
first coordinate of the tangent vector under the flow of the variational equation is first 
positive, then negative, then positive, and finally ends up negative because of lemma [5j 
These three sign changes mean that the vector has passed through the vertical at least 
three times, and hence we have an unstable standing wave. 

Remark 2. The solution shown in Fig. [6]^a) (see the corresponding phase-portrait in 
Fig. IA2( a)) belongs to the class of solutions discussed in Corollary HJ Nevertheless, 
the symmetric solution corresponding to point C in Fig. |5] (cf. Fig. EJc) and its phase- 
portrait [A2tc)) does not belong to the same class of solutions. In the phase-space, the 
orbit of the latter solution is sketched in Fig. [TiTa). This is similar to the configuration 
of the orbit in Corollary HI which is a symmetric (about the w-axis) solution in the phase 



space with the properties ^ < |, y < and = (—u ,y ), but here u < y^. 

In this case, the inner orbit lies outside the heteroclinic orbit of the inner system. In 
this configuration, we have that Q = 1 and P > 2. This implies that the above method 
cannot be used to prove the instability of such a solution. In Fig. [TITb). we plot the 
positive spectrum > of the operator D + for unstable symmetric solutions in Fig. 
[51 where one can see that P = 3 only in some interval. We conjecture that the break 
down of our method occurs at the points where the branch containing point B and that 
containing points D and E emerge with the main branch corresponding to symmetric 
solutions. 



4-3. Higher order excited states 

We remark that there are certain choices of V and u, and u such that the inner part 
of the superimposed phase portrait will be part of a periodic orbit of the inner system. 
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(a) 




0.65 



0.7 



0.75 
CO 



0.85 



(b) 

Figure 11. (a) A sketch of symmetric orbits with P > 2 and Q = 1. (b) The same as 
in Fig. |9l but for unstable symmetric solutions in Fig. [5] 



This applies for example if we keep the (uo,yo) and (ui,y\) = (— «o,J/o) as in corollary 
HI we can allow the trajectory to follow the periodic orbit for as many periods as we 
like, and can calculate the net effect on the stability of the standing wave as a whole. 
One can show that if the base point travels along the inner orbit described by (1161) from 
(ui, yi) for exactly one period, the solution to the variational equation will have exactly 
2 more zeros, as will the original solution. Thus for each period P = Q = 2, and so 
we still have an unstable excited solution. In this way we are capable of constructing 



unstable excited solutions with 2k + 1 zeros for an arbitrary k. Figure 10(b) is a phase 
portrait for such a standing wave solution. 

It remains to determine which choices of V, w, and Uq that will produce such an 
excited state. Using the geometry of the superposed phase portraits we have that such 
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an excited steady state solution will exist when 



< 



and u < 



V > \ \ an d u o < 



V-W 2 -2u 2 



V 
2 ' 



or 



(27) 



The above technique can be used to create a higher order excited state whenever the 
inner part of the phase curve is part of a periodic orbit. We can simply follow the period 
as many times as is necessary and the higher order state will have the same difference 
of P and Q. Thus if we begin with an unstable state where P = Pq and Q = Qq being 
the initial values for P and Q and the inner part of the orbit is periodic (for another 
example c.f. figure(8(a) )), we can produce higher order unstable states with P = 2k + P 
and Q = 2k + Qq, and for all integer values of k we will still have an unstable state. 



5. Asymmetric states 




(b) 



Figure 12. (a) An asymmetric orbit with P > 1 and Q 
with P = 2 and Q = 1. 



0. (b) An asymmetric orbit 



The techniques described and used in the section above cannot be applied to 
the asymmetric states considered in this paper because such states do not satisfy the 
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condition in theorem [TJ As particular examples, we consider two solutions with Q = 
and Q — 1. 

For positive solutions Q = 0, when y < |, asymmetric orbits will be present. An 



example of the phase portrait of a positive asymmetric solution is in Fig. 12(a) For 



this case, it can be calculated that P > 1. The profile shown in Fig. [3]^c) with its 
phase-portrait depicted in Fig. IAl( c) belongs to this class of asymmetric solutions. As 
we know numerically that asymmetric positive solutions are stable (see Fig. [3]), it is 
expected that Theorem [1] should not apply here. It is then interesting to know whether 
the stability can be established analytically. 

Numerically (not shown here) we observed that P = 1 for the asymmetric positive 
state. If one could show analytically that P = 1 indeed, one would be able to prove 
the stability using the theorem presented in [3Tj . Nonetheless, the solution stability is 
expected and predicted using the method of [19J. 

Excited asymmetric orbits where P = 2 and Q = 1 can occur for all configurations 
of V and w. Such an orbit is described by the following properties. If y < §, then for the 

v 



point of departure (w ,?/o), we have that uq < y ^ and yo < 0. The orbit then travels 
through the y axis in the phase plane and reconnects with the outer orbit at the point 
(ui, yi) where u\ < and y\ > 0. Moreover, in this configuration of V and u, the only 
type of inner orbit that can produce such an asymmetric solution happens to be periodic 
and symmetric about the line y = u in the phase plane. Thus (ui,yi) = {—Uq, —yo), and 
the original vector tangent to the homoclinic orbit (/3 ) gets mapped by the variational 
flow of the inner orbit to the tangent space of the homoclinic orbit at the point {u\, y\). 
Thus we know 

^ (^(P )( Uuyi ))( U ,y) f -k 
x™ |$ (^(^°)(ui,yi))(«,i/)| ^ kVV^ 

That is the limit under the variational flow of the vector which was originally tangent 
to the homoclinic orbit is the stable subspace of the variational flow at the critical point 
at the origin. We can therefore only conclude that P is exactly 2. If we are in the case 
where y > |, and < u < \J\, then we can exploit the same symmetry conditions 
as above for such asymmetric solutions, even concluding the limit of the variational 
solution under the tangent flow as the stable subspace of the critical point at the origin. 
An example of one such orbit is in figure 12(b) , which is a sketch of the profile shown 

in Figs. |6]^d) and (e). If however, y > | and uq > y^, then no such periodic orbits 
exist, i.e. there is no connected inner orbit where (u ,y ) and could have the 

aforementioned properties. 



6. Summary 

We have considered a nonlinear Schrddinger equation with a non-uniform nonlinearity 
coefficient. In particular, we have investigated a Schrddinger equation with self- 
defocusing nonlinearity bounded by self-focusing one. The present work extended the 
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results of [211 [13 US] • We have established analytically the instability of symmetric states 
beyond a critical norm through the application of a topological argument developed in 
[13] . Even though the technique does not definitively establish instability in the case 
of the asymmetric states considered in this paper, we have numerically established the 
stability of some positive asymmetric states. The analytical (in) stability of asymmetric 
higher-order modes remains an open problem, which is proposed to be studied in the 
future. 

Appendix A. Phase portraits 

In Fig. [31 the solution profiles corresponding to points A-D in Fig. [2] are shown in the 
physical space. The corresponding phase-portraits of the solutions are presented in Fig. 
IA11 The phase-portraits of solutions shown in Fig. [6] are depicted in Fig. IA2I 




0.2 0.4 0.6 0.8 1 '0 0.2 0.4 0.6 0.8 1 

u u 



(c) (d) 

Figure Al. Phase portraits of solutions at points indicated as A-D in Fig. [2] (see also 
Fig.©. 
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